Physics of interface: Strongly Correlated Barrier with chemical modulation 
sandwiched between two metallic planes. 
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Barrier planes described by the Ionic Hubbard model and sandwiched between metallic planes on 
both sides are studied using unrestricted Hartree Fock. For zero onsite correlation, the presence of 
the metallic interface generates an additional gap in the energy spectrum away from half filling, if 
the chemically modulated potential in the barrier planes exceed a critical value. There is reentrant 
behaviour and an insulator-metal-insulator transition as we tune onsite correlation for fixed site 
potential. The metal is thus able to penetrate the barrier planes due to proximity effect, in a system 
that is otherwise an insulator throughout. 

PACS numbers: 73.20.-r, 71.27.+a, 71.30.+h, 71.10.Fd, 73.21.b, 73.40.C 
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INTRODUCTION 

The emergence of unexpected new properties has 
made the physics of interface an active research area 
for some time now[JQ. The interface of a band in- 
sulator and Mott insulator was shown to have metallic 
properties which becomes superconducting on low- 
ering the temperature [H, Hj]. Theoretically the subject 
has been studied using different methods like restricted 
Hartree Fock@, two site DMFTQ, single site DMFT 
0, HI and Lancsoz @ . 

In this letter we have studied a metal-quasi 2D barrier- 
metal hetcrostructurc. The barrier planes are described 
by the Ionic Hubbard Model, which describes a band in- 
sulator with onsite correlation. Thus the barrier planes 
describe a band insulator in the absence of correlation U, 
and a Mott insulator in the absence of chemical modula- 
tion for U > U c , where U c depends on the dimensionality 
of the system. The ionic Hubbard model has been widely 
studied by variousgroups using various techniques like 
single site DMFT [jj quantum monte carlo[ll| and 
cluster DMFT Q3 in 2D, DMRG in ID [HI and two site 
DMFT in quasi 2D At half filling a single site DMFT 
solution which assumes a paramagnetic solution yields a 
metallic phase which separates the Mott insulator and 
band insulator phases 0, A solution which incorpo- 
rates long range anti-ferromagnetic ordering in quasi 2D, 
using two site DMFT, however shows that the system is 
insulating for all interaction strengths. 

The Mott insulator-band insulator heterostructure has 
been studied using restricted Hartree Fock inhomo- 
geneous two site DMFT Q and in quasi one dimen- 
sional lattices using Lancsoz method @. The metal- 
Mott insulator-metal heterostructure has been studied 
@, B EH while the role of disorder on the Mott planes 
has also been investigated by us [l6[ . 



Spatial variations in the planes and along the z di- 
rection play a crucial role in determining the correct 
ground state of a highly inhomogeneous systems like het- 
erostructures. The in plane variation is totally ignored 
in a one site inhomogeneous DMFT(IDMFT)[lJ-l2(i| ap- 
proach. The two site IDMFT method though capable of 
capturing planar long range order like planar antifcrro- 
magnetism, is insufficient to capture the regime where 
the two orders start to cancel each other leading to short 
range ordering. With this in mind we have adopted the 
method of unrestricted Hartree Fock. The Hartree Fock 
method has been used successfully to describe the ground 
state of the ionic Hubbard model even in ID and was 
shown to be in excellent agreement with real space RG 
calculation 21 1 . 



Model and Method 

The Hamiltonian for the system is 

7~t = — tij c lacr c jao- ~ ^ [ c \ao c ia+ltj + h.C.] 
ijaa iaa 

iaf T" ^zaj, ) 



(1) 



Here the label a indexes the planes, and the label i in- 
dexes sites of the two-dimensional square lattice in each 
plane. The operator c\ aa {c iaa ) creates (destroys) an 
electron of spin a at site i on the plane a. We set the in- 
plane hopping t" to be nearest neighbor only, and equal 
to t, the hopping between planes, so that the lattice struc- 
ture is that of a simple cubic lattice. We take W a = W 
for the B sites and W a = for the A sites. We also take 
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FIG. 1: Energy spectrum for U — for a) W — 1, b) W = 
4, c) W = 10, d)W = -4, e)W = -10. f) Evolution of the 
multigaps with increasing W for [7 = 0. 



U a = U for the barrier planes, and zero for the metallic 
planes. The chemical potential fx is calculated by tak- 
ing the average of the N/2 th and the N/2 + 1 th energy 
level. The plane index a = 1, m correspond to the metal- 
lic layers. All the other a values in between correspond 
to the barrier planes. For the quasi 2D geometry taken 
by us total number of sites N = m x L 2 , where m is the 
total number of layers and L 2 is the number of sites in 
the planes. 



Calculated quantities 

We have calculated the energy spectrum, charge and 
spin profile, and the optical conductivity. The charge 
at site i is given by: Cj = n^-f + H,,^. The spin at a 
particular site is given by Si = \(rii^ ~ n^)). The op- 
tical conductivity is calculated using the Kubo formula. 
In this paper, a zz (ui) is calculated in units of ire 2 /haQ. 
For details of the calculation procedure see our previous 
work 22 , [23| . We have performed finite size scaling of 
each of our transport results. The value of a zz {ui) is cal- 
culated for multiples of a small u — tu r , where uj r is twice 
the lowest uj that we can access for a particular system 
size. 



Analysis of our results 

Fig la,b,c,d,e and f shows that different types of 
sites emerge as we systematically tune W. W posi- 
tive/negative corresponds to repulsive/attractive poten- 
tial/wells in the barrier planes. In the absence of U any 
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FIG. 2: a. Gap at half filling for W — 1 with increasing U. 
There is a prolonged gapless phase that extends from U — 1 
to U — 7. b. Gap at and away from half filling for both 
spins at U = 7 with W. c. Energy of the two spins plotted 
separately along with the total energy, d. Charge profile at 
W = 1 and U = 3. 



non-zero W leads to a charge density wave(CDW) profile, 
with a symmetry in the gap about the origin. Increas- 
ing W beyond a certain threshold value, an additional 
gap opens up away from half filling. The value of the 
gap which would have been identical to W in the ab- 
sence of the metallic planes gets reduced in the presence 
of the metallic planes. This gap is due to the fact that 
the barrier planes induce a density modulation even in 
the metallic planes due to proximity effect. Thus the 
entire system has a rather inhomogeneous alternate den- 
sity modulation. As we increase W further at around 
W g = 4, a second gap opens up. For general m, the fill- 
ing at which the second gap opens up is given by 0.5 + 
sign(W)*(l/m). Above this value W g , the lowly occupied 
sites in the barrier planes (B sites) become very different 
from the lowly occupied sites in the metallic planes and a 
new band emerges. For negative W the second gap opens 
at W = -4. The states below this new band gap now cor- 
respond to the B sites from the barrier planes and the 
states above the gap correspond to lowly occupied sites 
from the metallic planes. Due to -ve value of W the sys- 
tem can now lower energy by doubly occupying the B 
sites. The gap at half filling saturates to a fixed value 
as we keep on increasing W. This is because the higher 
energy B sites (both in the barrier and those induced in 
the metallic planes) are now more or less avoided by the 
electrons. 

In Fig 2a, we show the modulation of the gap at half 
filling for W = 1 as we increase U. The band gap is sup- 
pressed by U down to zero, after which there is a large 
gapless region. The low U region corresponds to a CDW 
phase, where U has started to compete with the chem- 
ical modulation with the latter being dominant. The 
large gapless region corresponds to a highly inhomoge- 
neous phase, with no spin order till U — 6 and with 
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FIG. 3: a. Gap at and away from half filling for W = 4 with FIG. 4: Optical conductivity results ct(oj) vs ui for a. {7 = 1, 

increasing U. b,c,d. Charge profiles for W — 4 and U — 2, 3, 6 b. (7 = 3, c. U = 6, d. U = 7. 

respectively. 



gradually diminishing but non zero charge order. Above 
[7 = 6, spin order starts to pick up gradually, though 
there is still a remnant weak charge order also. 

Fig 2b shows the variation of the two gaps in the energy 
spectrum for U = 7, for both spins with variation in 
W . In the absence of W , the opening up of two gaps in 
the e nerg y spectrum was first shown by us in a previous 



work [15j. For positive values of W, till W = 1, both 



the gaps for both spins fall monotonically with increasing 
W . There is small spin asymmetry which is a signature of 
spontaneous symmetry breaking of the up and down spin 
sectors for such a system. This behaviour has already 
been observed and reported by Cabib and Callen [24[ for 
an extended Hubbard (t — U — V) model. 

While for positive W the B sites are avoided by the 
electrons, The repulsion of the electrons to the metal- 
lic planes takes place uniformly for both up and down 
spins because there is very small charge modulation in 
the metallic planes and since there is U on all the barrier 
plane sites, they would all like to push electrons away to 
the metallic planes, without paying the cost of U . 

Fig 2c shows the plots of energy for both spin species 
and the total energy as we modulate W for U = 7. The 
total energy increases monotonically with increasing W. 
The energy for the two spin species are different for the 
entire parameter range where there is both charge and 
spin order. This is a signature of the spin asymmetry 
reported earlier |24j. 

2d shows the plots of the charge profile for [7 = 7. The 
overall symmetry of the profiles about the central plane 
and the explicit z dependence are seen. In this regime 
where W has just killed the gap due to anti ferromagnetic 
long range order due to [7, there is already a strong plane 
dependent CDW ordering. 

In Fig 3a we have shown the variation of the two gaps 
that opens up for W = 4 with increase in [7. We see that 
the two gaps behave in very contrasting manner with 
increasing U. While the gap at half filling falls off mono- 



tonically to zero, which basically implies that the charge 
ordering induced in the metallic planes by W is killed by 
[7. However the second gap for positive W > 4 shows a 
non monotonic variation with [7. It first increases with 
increasing U and then after a certain value of U starts to 
decrease. This is because for low values of [7, the elec- 
trons in the A sites of the interface layer start to deplete 
by moving to the metallic planes. The occupancy of the 
B sites however is relatively unaffected in this regime. 
Thus the second gap which senses the difference between 
the topmost level of the second band and the lowest level 
in the topmost band increases. However beyond a certain 
threshold, the B sites in the barrier planes start to get 
occupied and then the gap starts falling. 

Figs 3b,c,d show the charge profile for U = 2,3,6 re- 
spectively for W = 4. While for U = 2 and 3, the gap 
at half filling decreases but is still non zero, for [7 = 6 
the gap at half filling has already become zero. This we 
can see clearly from the figures, where for [7 = 2 and 3, 
there is inhomogeneous CDW ordering in all the planes. 
For U = 6 however the CDW ordering is destroyed in the 
metallic planes. 

Figs 4a,b,c,d show the plots of optical conductivity 
versus uj for [7 = 1,3,6 and 7 respectively for W = 1. 
As we tune U keeping W fixed, the system exhibits 
reentrant behaviour in the form of an insulator-metal- 
insulator transition. As the gap at half filling gets killed 
initially around U = 1, a residual slightly weakened CDW 
ordering persists in the system all the way upto [7 = 6. 
The A sites in the barrier start to deplete and the B sites 
are relatively unaffected for low [7. As we increase U fur- 
ther, the system becomes a metal as seen in U = 3 and 
6. The metal at U = 3 is a much weaker one compared 
to the one at U = 6, as the charge ordering is fully de- 
stroyed for [7 = 6. The weak metallic behaviour for U = 
3 inspite of the presence of charge ordering in the system 
is a consequence of weak tunneling of electrons between 
the sites. For [7 = 7, a small inhomogeneous spin order 
develops in the system, which drives the system to an 
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insulating phase. 



CONCLUSION 



A barrier described by the ionic Hubbard model sand- 
wiched between two metallic planes has been studied us- 
ing unrestricted Hartrce Fock. Increasing the strength of 
the site potential modulation, beyond a critical strength 
W c opens a new gap away from half filling at a value that 
depends on the sign of the modulation. Below W c onset 
of correlation reduces the gap to zero gradually. Above 
W c while the gap at half filling follows the same pattern, 
the new gap shows non monotonic behaviour with in- 
creasing U. The system shows reentrant insulator-metal- 
insulator transition with increasing U for fixed W = 1. 
Spin asymmetric behaviour is noticed when both charge 
and spin order is present in the system. 
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